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Branching within branching I: The extinction 
problem 


Gerold Alsmeyer and Soren Grottrup 


Abstract We consider a discrete-time host-parasite model for a population of cells 
which are colonized by proliferating parasites. The cell population grows like an or¬ 
dinary Galton-Watson process, but in reflection of real biological settings the mul¬ 
tiplication mechanisms of cells and parasites are allowed to obey some dependence 
structure. More precisely, the number of offspring produced by a mother cell de¬ 
termines the reproduction law of a parasite living in this cell and also the way the 
parasite offspring is shared into the daughter cells. In this article, we provide a for¬ 
mal introduction of this branching-within-branching model and then focus on the 
property of parasite extinction. We establish equivalent conditions for almost sure 
extinction of parasites, and find a strong relation of this event to the behavior of par¬ 
asite multiplication along a randomly chosen cell line through the cell tree, which 
forms a branching process in random environment. In a second paper G), the case 
when parasites survive is studied by proving limit results. 
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1 Introduction 


The discrete-time branching-within-branching process (BwBP) studied in this paper 
describes the evolution of generations of a population of cells containing proliferat- 
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ing parasites. In an informal way, its reproduction mechanism may be described as 
follows: 

(1) At time n = 0 there is just one cell containing one parasite. 

(2) Cells and their hosted parasites within one generation form independent re¬ 
production units which behave independently and in the same manner. 

(3) Any cell splits into a random number N, say, of daughter cells in accordance 
with a probability distribution (pk)k> o- 

(4) Then, given N, the hosted parasites, independently and in accordance with 
the same distribution, produce random numbers of offspring which are then 
shared into the daughter cells. 

(5) All cells and parasites obtained from a cell and its parasites in generation n 
belong to generation n + 1. 

We are thus dealing with a hierarchical model of two subpopulations, viz. cells and 
parasites, with an entangled reproduction mechanism. The hierarchy stems from the 
fact that cells can survive without parasites but not vice versa. 

Proceeding with a more formal introduction, let V denote the infinite Ulam- 
Harris tree with root 0 and N v the number of daughter cells of cell v £ Y. The 
(A v )v 6V are independent and identically distributed (iid) copies of the No-valued 
random variable N with distribution (pk)k>o and finite mean v, viz. P (N = k) = pk 
for all k £ No and 

v = E N < °°. 

The cell population thus forms a standard Galton-Watson tree (GWT) T = |J,, 6No ^ 
with To = {0} and 

T „ := {vi ...v„ £ V|v!... v„_! £ T„_i and 1 < v„ <A Vl ... Vn _j} 

(using the common tree notation vi ...v„ for (vi..... v„)j. Consequently, defining 

X ■■= #T„ = £ Av (1) 

vgT„_| 

as the number of cells in the n th generation for n £ No, the sequence (%)„>o forms 
a standard Galton-Watson process (GWP) with reproduction law (pk)k >0 a 'id repro¬ 
duction mean v. For basic information on Galton-Watson processes see HU El. 

Let Z v denote the number of parasites in cell v £ V and T* the set of contaminated 
cells in generation n £ No with cardinal number PG , so 

T* := {v £ T„ : Z v > 0} and Sff := #T*. (2) 

We define the number of parasites process by 

■■= E Z v> » e No. 

veT„ 

After these settings, the BwBP is defined as the pair 
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(T„, (•Z v )vgT„) )1 >0 


and now clearly seen as a description of the generations of any given population of 
cells and the number of parasites hosted by them. 

As informally stated above, parasites at different cells are assumed to multiply 
independently of each other, whereas reproduction of parasites living in the same 
cell v is conditionally iid given the number of daughter cells of v. More precisely, 
let for each k £ N 


y (•>*) 

A i,v 


:= X, 


(U) 


.A, 


(k,k) 


i £ N, v G Y, 


be iid copies of an Ng-valued random vector := (x ,] ’ k \... . The fam- 

ihes k £ N, are assumed to be mutually independent and also inde¬ 

pendent of (Av) v gv- They provide the numerical description of the reproduction and 
sharing of the parasites living in the cell tree. In detail, if the cell v £ V has k £ N 

(j ]A 

daughter cells vl, ...,vk, then A A’ , 1 < j < k, gives the number of progeny of the 

i th parasite in cell v which go in daughter cell v j. The sum over all entries in xj*'^ 
gives the total offspring number of this parasite. The number of parasites in the cells 
are thus recursively defined by Z 0 := 1 and 


Zvf := E %„=*} E x if = E X iy Ny) - j£ 

k>j i=l i=l 


(3) 


(j 

where, by convention, ’ : = 0 if j > k. We further define 


Hk ■= E 


for j,k e N and 

7 := = E Pk E Hk 

k> 1 y=l 

as the mean number of offspring per parasite, which is assumed to be positive and 
finite, i.e. 

0 < y < °o. (Al) 

This naturally implies jl i^ < °° for all j < k, and P (N = 0) < 1. To avoid trivial 
cases, we further assume that 

pi = P (N = 1) < 1 and P( = 1) < 1. (A2) 

If the first assumption fails, the cell tree is just a cell line and (J'X n )n>o a standard 
GWP with reproduction law ,.X'(X nA l ), while failure of the second assumption en¬ 
tails the number of parasites in each generation to be the same, thus ST* = a.s. 
for all n £ No, or S'* = So eventually. To rule out the simple case that every con- 
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Fig. 1 A typical realization of the first three generations of a BwBP. 


taminated daughter cell contains only one parasite (as the root cell 0), we further 
assume that 


Pk > 2) > 0 for at least one (j,k), 1 < j < k. (A3) 

For each generation n = 0,1,2,..., let 


Tn •— ^n, 1 1 «^n,2) • • ■) 


denote the infinite vector of cell counts with a specific number of parasites, so 
gives the number of non-infected cells and for k > 1 the number of cells with 
exactly k parasites in generation n. Then (77i)neNn i s a multi-type branching process 
(MTBP) with countably many types. The individuals are the cells and the type of 
each cell is given by the number of parasites it contains. 

In situations where our BwBP initially has sk cells of type k for k G I Io and thus 

a total number . () ,v, of cells, we write P s for the underlying probability measure, 

where 

s = (s 0 ,si,...) e N := {(*,-),>0 G Nq|x, > 0 finitely often}. (4) 

Further, we write P Z) z G No if the process starts with one cell and z parasites, i.e. 

F z (% = l,Z 0 =z) = l. 

The corresponding expectations are denotes as usual by E s and E-, respectively. 
Indices are dropped in the standard case, viz. P = Pi and E = Ei. 
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MTBP’s with finite type-space are well-studied with results transferred from the 
classical theory of GWP’s, see no Chapter V], (29] Chapter 4] or the monography 
by Mode ll37l . If, on the other hand, the state space is infinite (countable or uncount¬ 
able), a variety of behaviors may occur depending on the reproduction mechanism 
of types. For example, the branching random walk, an object of long standing inter¬ 
est HU 0 El M US M EU Ell IH , is obtained when type-space transitions are of 
independent additive kind. MTBP’s are also studied in J8] [22 ] [3T] [36] [38] 1391 14011 . 
and we refer to Kimmel and Axelrod 1)33] Chapter 7] for a series of examples with 
applications in biology. 

Various examples of branching-within-branching models have already been stud¬ 
ied in the literature, the first one presumably by Kimmel l32l who considered binary 
cell division in continuous time with symmetric sharing of parasites into the daugh¬ 
ter cells. A discrete-time version of his model, and in fact a special case of ours, with 
possible asymmetric sharing of parasites was studied by Bansaye El and later ex¬ 
tended in ITU by adding immigration of parasites and random environments. The 
latter means that parasites in a cell reproduce under the same but randomly cho¬ 
sen distribution. Subsequently, efforts have been made to generalize the underlying 
binary cell tree of the afore-mentioned models to arbitrary Gabon-Watson trees in 
both, discrete and continuous time. The greatest progress in this direction has been 
achieved by Delmas and Marsalle in ll24l and together with Bansaye and Tran in 
El . Both articles consider a random cell splitting mechanism and asymmetric shar¬ 
ing but make ergodic hypotheses that rule out the possibility of parasite extinction 
which is the focus in the present article. Let us further mention work by Guyon |[26l 
on another discrete-time model with asymmetric sharing and by Bansaye and Tran 
El on a bifurcating cell-division model in continuous time with parasite evolution 
following a Feller diffusion. There, the cell division rates depend on the quantities of 
parasites inside the cells and asymmetric sharing of parasites into the two daughter 
cells is assumed. 

Beside the model in |T2l . our model also comprises the one of type-A cells stud¬ 
ied in 0, where these cells produce daughter cells of either the same type A or of 
type B, and the sharing mechanism of parasites in a type-A cell may depend on the 
number of type-A and type-B daughter cells. As another special case, we mention 
the situation of multinomial repartition of parasites. Here all parasites multiply in¬ 
dependently in accordance with the same offspring distribution so that the number 
of parasites process {3f n ) n >(t forms a standard GWP. After parasite reproduction, 
a cell divides into a random number N of descendants, and the offspring of each 
of its hosted parasites chooses independently the i lh daughter cell with probability 
qi(k) £ [0, 1] if N = k. Thus, 

£ xM L x(M) (5) 

7=1 

for all k £ N, and given Y}j= \ X^' k > = x, the vector .. ,X^ k ^) has a multi¬ 

nomial distribution with parameters x and q\(k),... .qiJk) £ [0,1]. Replacing ([5]> 
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with 

x U, k ) £ x (L1 ) for all 1 < j < k < °°, 

we arrive at a case in which the number of offspring shared by a parasite into a 
daughter cell is always the same and thus independent of the number of such cells 
produced. This implies that along an infinite cell line the number of hosted para¬ 
sites forms an ordinary GWP. As yet another specialization, one may consider the 
situation when 

xU, k ) = o a.s. for all 2 < j < k < 

while are independent and positive-valued. Hence, starting with a 

single cell, parasites are only located in the leftmost cell 1„ = 1...1 (n times) in each 
generation n and (3f n ) n > o forms a GWP in random environment (GWPRE), with 
the number of daughter cells of each 1„ forming the (iid) random environment and 
hence determining the offspring distribution in each generation (see cni eh for the 
definition of a GWPRE). 

Regarding our model assumptions, let us note that it is not for pure mathemati¬ 
cal generality when allowing cell division into more than two daughter cells. In the 
standard case, one of the two daughter cells after cell division may be viewed as 
the original mother cell which accumulates age-related damage throughout its repli¬ 
cation phases and eventually loses the ability for cellular mitosis so that cell death 
occurs. This phenomenon, called cellular senescence, has been discovered recently 
even for several single-celled organisms (monads), see Stephens 11421 . Genealogical 
aspects may be studied with the help of Galton-Watson trees when counting all cells 
stemming from a single cell during its lifetime and interpreting them as the suc¬ 
ceeding generation. Since the infection level of the mother cell changes during its 
lifetime, different numbers of parasites in the daughter cells are to be expected, thus 
justifying the assumption of asymmetric sharing of parasites. As another reason for 
this assumption is the fundamental biological mechanism to generate cell diversity, 
see Jan and Jan f30l and Hawkins and Garriga |27| . For example, a stem cell uses 
asymmetric sharing to give rise to a copy of itself and a second daughter cell which 
is coded to differentiate into cells with a particular functionality in the organism. 

The purpose of this article is to establish equivalent conditions for almost sure 
extinction of parasites. The limiting behavior of the BwBP will be studied in a com¬ 
panion paper 0. 


2 The associated branching process in random environment 

The first step towards an analysis of the BwBP is to identify a certain infinite ran¬ 
dom cell line through the cell tree T and to study its properties. This approach was 
first used by Bansaye in m who simply picked a random path in the infinite bi¬ 
nary Ulam-Harris tree representing the cell population. Since the cell tree is here a 
general GWT and thus random, we must proceed in a different manner already in- 
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Fig. 2 Typical realization of a spine in the size-biased cell tree. A spinal cell is shown as a red 
square □, all other cells as a circle Q- 


traduced in (6). In fact, we pick the random path (V„)„>o in V according to a certain 
size-biased distribution instead of uniformly. 

To give details, let (C n ,T n ) n >o be a sequence of iid random vectors independent 

of (/V v )vev and (X^*’* )*>i j>i veV- The law of T n equals the size-biasing of the law 
of N, i.e. 


P (T n = k) 


kpk 

v 


for each n £ No and k £ N, and 


P(C„ = l\T n =k) = y 
k 

for 1 < I < k, which means that C„ has a uniform distribution on {1,..., k} given 
T n = k. Now, (V„)„>o is recursively defined by Vo = 0 and 


Yi •— Vn— \C n — i 


for n > 1. Then 

0 =: Vo —V Vi —)■ V 2 —► • • • —y V n —>... 

provides us with a random cell line in V (not picked uniformly) as depicted in Figure 
[2] The resulting path can be viewed as a so-called spine in a size-biased tree, each 
spine cell representing a “typical” cell of its generation in the ordinary cell tree. The 
concept of size-biasing in the branching context goes back to Lyons et al. in l35l . 
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who used it to give alternative proofs of classical limit theorems for GWP’s, and we 
refer to their work for a detailed construction of a spinal GWT. 

Concentrating on the number of parasites along ( V n ) n >o, we have Zy 0 = Z 0 and, 
recursively. 


z v„+i = 


“ k Zv " t ■ 

E E 1 {C„=j, T„=k} E X i.v„ = 

*=u=l i= 1 


Z Vn 

L x < 

i= 1 


(CnJn) 

UV„ 


for n > 0. Hence, given ( C n . T n ), all parasites in generation n produce offspring 
independently with the same distribution Since the ( C n ,T n ) are iid and 

independent of (^/*’^)t>t i>i,veV> the process of parasites along the spine (Zy n ) n >o 
forms a branching process in iid random environment (BPRE) as defined e.g. in i TOl 
ED. We summarize this observation in the following lemma (see also |6] Subsection 
2 . 1 ]). 

Lemma 2.1 Let (Z' n ) n > o be a BPRE with Z 0 ancestors and iid environmental se¬ 
quence A \= (A n ) n >o taking values in {jSf(X^"’^)|l < j < k < °°} and such that 

v(Ao = A?(xW)) 

for all 1 < j < k < °°. Then (Zy n )„>o and (Z' n )„>o are equal in law. 

Proof. It suffices to point out that, by assumption, Z' 0 = Z 0 a.s. and 

P (Ao = f = \ ■ ^ = P (Co = j,To = k) 

for all 1 < j < k < °°, i.e., both processes pick their reproduction law in each gener¬ 
ation by the same random mechanism. □ 

The BPRE (Z')„>o with environmental sequence A is called hereafter the asso¬ 
ciated branching process in random environment (ABPRE). It is one of the major 
tools used in the study of the BwBP, and the following proposition provides a key 
relation between this process and its ABPRE. 

Proposition 2.2 For all n.k.z £ No, 

P z (Z' n =k) = V ”E z 3T n k and P, (z' > 0) = V'” E Z £T. (6) 

Proof. By definition of the involved random variables, we find that 


-^zC^vlO ^0? • • • 1 ^n— 1) 


n— 1 


n pl 

i =0 


V”P z (V n 


v, 7o = ko,...,T n -i =k„-\) 


and 
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E z( Z v — ^|-^v |0 — k 0 , • ■ • — k n — l) 

= P-( z v = k\v n = v, T 0 = k o,..., 1 = 

for all vertices v = Vi... v„ and ko > vi,... ,k„-i > v„, where v |0 := 0 and 

v| / := V| ... v ; for 1 </<;; = |v|. Hence it follows by summation over the Ay that 

P z (Z v =k) = v" P;(Z V = k,V„ = v) 

for all v £ V with |v| = n, which in turn leads to 

v" P : (Z y „ = k) = £ v" P,(Z V =k,V n = v) 

|v|=« 

= X E -( Z v=£) = X E Z 1L {Z,=H = E Z-4 

|v|=n | v|=« 

and proves the first assertion in view of Lemma |2T| The second equation then fol¬ 
lows by summation over all k > 0 . 

We finish this section with a quick review of some relevant facts about the BPRE, 
relevant references being 1121 [3] [4] O[TO][14] [23] [25] [TT] |43] IS] 1451 . For n C N and 

s e [ 0 , 1 ], 

E(s z " |A) = gAo°-OgA n _,(s) 

is the quenched generating function of Z! n with iid g\ n and gx defined by 
gx(s) := E (^‘ |A) = A)= X^« s ” 

n> 0 

for any distribution A = (A„)„>o on No- Moreover, 

£ 4 ,( 1 ) = ez; = x v ExU ' k) = = l < (7) 

1 <j<k V 

where y= EiZ’j. It is also well-known that (Z')„>o survives with positive probability 
iff 


E1 °g t ?A 0 ( 1 ) > 0 and Elo S (i-gAoW) <°°> (8) 

see e.g. EBBD, and recall that y < °° is assumed by ( |A 1 [ . Furthermore, by ( |A3[ ), 
there exists 1 < j < k < °° such that pk > 0 and P(X^'^ ^ 1) > 0, which ensures 
that Ao with positive probability. The ABPRE is called supercritical, critical 
or subcritical if Elogg( io (l) >0, = 0 or < 0, respectively. The subcritical case 
further divides into the three subregimes when Eg( Vo (l)logg( io (l) < 0,= 0, or > 0, 
respectively, called strongly, intermediate and weakly subcritical case. The quite 
different behavior of the process in each of the three subregimes is shown by the 
limit results derived in ll25l . 
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Focussing on the subcritical case hereafter, thus Elogg^ 0 (l) < 0, we point out 
the following useful facts. If Eg^ o (1) log(1) < 0, then the convexity of 0 i —> 
Eg^ o (l) e implies that 

E ^o(!) = 

Under the assumptions 


P(Z',<C) = 1 and P (0 < g^(l) < e) = 0 (AsGe) 

for suitable constants C > 0 and e > 0, Geiger et al. ]25l Thms. 1.1-1.3] showed 
that 

P(Z' > 0) ~ cn~ K ^ inf ^Eg^ o (l) e ^ as n —* °° (9) 

for some c £ (0,°°), where 


! 0 if (1) log g' A ^ (1) < 0 (strongly subcritical case), 

2 if (1) log g' An (1) = 0 (intermediate subcritical case), 

2 if Eg^ o (1) log(1 ) > 0 (weakly subcritical case). 

The condition |AsGe]l can be strongly relaxed for the asymptotic relation of the 
survival probability in (|9]i to hold, see for example [25; 46j, but is enough for our 
purposes. 

A combination of ([6]) and ([9]) provides us with the asymptotic relation 


E3T = QsfjEskU)') 

in particular (with ( |AsGe[ > still in force) 

inf Eg{. (l) e < — if supEi^* < 
o<e<t ^ v „>i 


as n —> oo 


( 10 ) 


(ID 


3 The extinction problem 


The trivial fact that non-contaminated cells are unable to produce infected daughter 
cells shows that the process 


V ■= 

also forms a MTBP with type space {1,2,3,...}. As a consequence, the following 
extinction-explosion principle is easily obtained for our model. 
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Theorem 3.1 (Extinction-explosion principle) The parasite population of a BwBP 
either dies out or explodes, i.e. for all s G N 


— > 0) + P s (-S2« —> °°) — 1. 


Proof. Let s £ N with ,v, > 0 for at least one i > 1. We prove that s is a transient state 
for the Markov chain (7~ t *) n >o- Consider three cases: If P(/V = 0) > 1, then 

P ( 7 ;* ^ s for all n > 1 \T 0 * = s) > F(N = 0) EfSi > 0. 

Otherwise, if P(iV = 0) = 0 but P( Jfj = 0) > 0, then there exists a k £ N such that 
PkF(Lj =l X (j ' k) = 0) > 0 and thus 



Finally, consider the case when P(N = 0) = 0 and P= 0) = 0. Recalling ( |A2[ ), 
we then have P(i2fj > 1) > 0 and thus /JjtP(Ly=i > 1) > 0 for some k £ N. 
But this implies 



P (77 ^ s for all n > 1|77 = s) > p k P >1 > 0, 


that is, s is a transient state of the Markov chain (' Tf) n >t)■ Consequently, we infer 
that almost surely 

52 -%a —> 0 , 52 .% ,■ —> oo or, a fortiori, :7 nA —> °° for at least one i > 1 

i> 1 r>l 

as n —> oo. In any case, the parasite population a.s. dies out or tends to infinity. □ 
We denote by 


Ext :={£?„—> 0} and Surv := Ext c = {3f„ —► 


the event of extinction and of survival of parasites, respectively, and put 
Pj! := P s (-|Surv) and E* := E s (-|Surv), 

for s G N. Also, let P* and EJ for z € N have the obvious meaning. 

Having shown that (S n ) n ^o satsifies the extinction-explosion dichotomy, we will 
now turn to the process of contaminated cells and prove that, ruling out a degenerate 
case, survival of parasites always goes along with the number of contaminated cells 
tending to infinity. In other words, tends to infinity as n oo if this holds true 
for S„. The degenerate case occurs if all parasites sitting in the same cell send their 
offspring to the same daughter cell which formally means that, for each k > 1, 
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X^ k) = (0,...,0,X (; ' W ’ fc) ,0,...,0) 

for some unique j(k) £ {1, or, equivalently, 

p 2 (^r > 2) = o. 

Note that it is enough to consider a single root cell due to the branching property. 
Theorem 3.2 Let P(Surv) > 0. 

(a) If P 2 (^* > 2) > 0, then P* (Iff —> °°) = 1 and thus Ext = {sup ;l>0 If,* < °°j 
P z -a.s.for all z £ N. 

(b) If P 2 (^* > 2) = 0, then P* (iff = l fa. n > 0) = 1/or all z £ N. 

Proof Fix any zeN and consider first the easier case ( b). Note that P 2 (^* > 2) = 0 
implies P z ( Iff < 1 V n > 0) = 1. But since Surv = {Iff > 1 f.a. n > 0} P z -a.s., ( b) 
is proved. 

For the proof of (a), we use the Markov chain (Tf) n >o to show that (lff )„>o 
visits each t > 1 only finitely often, that is 

P z (l < ff < t infinitely often) = 0 (12) 

for all t > 1, hence P z (Um„_,.oo Iff = 0 or °°) = 1. But since Ext = { Iff -A 0} P,- 
a.s., (a) follows. 

Left with the proof of (p~2]), we define (with N given by (|4]>) 


Ak 


seN 


Y j s i = k, siylk 

i> 1 


c N 


for k > 1 and observe that, for n > 0, 


{If;=k} = {V&A,}U{Tf = (*,0,0,...)} P,-a.s. 


Since (*,0,0,...) £ N is a transient state by Theorem 3.1 we obtain 


P z (Iff = k infinitely often) = P z (7 ~f £ A^ infinitely often). 


Therefore, it remains to prove that the Markov chain (' Tf) n >o visits the set only 
finitely often with probability 1. For s £ A^, i.e. }_,■ s, = k and Sj > 1 for some j > 1, 
we infer with the help of the branching property 


P s (7;* i A k for all n> 1) 

> P s (ff > k for all n > 1) 

> P j(ff > 2) P(Surv) 2 Y[ (Pi(^t* > l)P(Surv)) 1 ' 

i:sj> 0 

> P 2 (^i* >2)P(t7/> lfP(Surv) i+2 
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and positivity of the last expression is guaranteed by our assumptions. Notice also 
that this expression is independent of the choice of s. Next, let To := 0 and 

T„ + i := inf [m > z n \T* £ A k } 

for n > 0 be the successive return times of (' T,*) n >o to the set A k , where as usual 
inf® := oo and T„+i = °° if T„ = °°. Then the previous inequality and the strong 
Markov property of {T~,*) n >o imply the existence of a constant c < 1 such that for 
all s £ A k and n > 0 

P Z (T„ + 1 -T„<°°| 7 ? n =S,T„ <°o) = P s (Ti < oo) < C < 1. 

Therefore, we infer upon iteration that 

P Z (T„ <°°) 

= ^ = S, 'In — ^n— 1 ‘n °°, ^n— 1 ^ 00 ) 

seAj. 

= ^ Pj(t w — T«_i < °°|— s,t h _i < °°)Pz(’7 "^_ 1 =S,T n _1 < 00 ) 
seA k 

< cP z (t„_i < <») < ... < c" _ 1 P ; (ti <°°) < c"” 1 

for each n > 1 and thereupon 

P- {T* £ A k infinitely often) = P, (t„ < °° for all n > 1) 

= F z (flK<“}) 

= limP-(T„<°°) < lime "” 1 = 0 

n —>oo ~ n— 

which completes the proof. □ 

So we have verified an extinction-explosion dichotomy for both, the process of 
contaminated cells and of parasites, and will now proceed with the main result in this 
paper which provides equivalent conditions for almost sure extinction of parasites. 
The proof will make use of the ABPRE introduced in Section [2] 

Theorem 3.3 

(a) If ¥ 2(^1 > 2) = 0, then P(Ext) = 1 if and only if 

ElogE(iFi|N 0 ) < 0 or Elog”P(iFi > O|jV 0 ) = °°. 

(b) If > 2) > 0, then the following statements are equivalent: 

(i) P(Ext) = 1. 

(ii) Efff < 1 for all n £ No- 

(iii) sup„ eN() E,9f < 00 . 
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(iv) V < 1 , or 


v>l, Elog^(l)<0 and inf E^ o (l) 0 < 


Proof, (a) If P 2 (^ > 2) = 0, then for all k > 1 with /?£ > 0 ther e exists at most one 


1 < /' < & such that Ppf- 7 -^ > 0) > 0, see before Theorem 3.2 As a consequence. 


{%i)n>o is 4 branching process in random environment, the latter given by the iid 
numbers of daughter cells produced by the unique cells which contain the parasites 
(their common law being the law of A 0 ). Hence, (f'Z n ) n> o dies out almost surely if, 
and only if, E(logE(,J°i |A 0 )) < 0 or Elog _ P(.2i > O|A 0 ) = °o (see e.g. fiTl ). 


(b) Suppose now Pi (> 2) > 0. 

“( 2 ) => {ii)” (by contraposition) Fix m £ N such that E (^*) > 1 and consider a 
supercritical GWP (S „)„>0 with .S’o = 1 and offspring distribution 


P (Si = k) = P(£T* =k), k£ N 0 . 


Obviously, 


for all k, n £ No, hence 


P (Sn > k) < P(^„; > k) 


limP(^;>0) 

n—¥ 00 


> lim P(S„ > 0) > 0, 

n —>00 


i.e. parasites survive with positive probability. 
“( 22 ) => ( 222 )” is trivial. 

“(Hi) => ( 2 )” Recall that liminf,,-*,*, STf = °° 
other hand, sup„ >0 E,5^* < °° implies 


a.s. on Surv by Theorem 3.2 On the 


00 > liminfE^* > E fliminf 

n —\ n —>00 / 


by Fatou’s lemma so that P(Surv) = 0. 

“(iv) => ( 2 ), ( 22 )” If V < 1, then E^* < E% = V n < 1 for all 22 G N. So let us 
consider the situation when 


v > 1) Elogg^(l)<0 and inf E^ (l) e < ^ 


is valid. It is here where the ABPRE (Z' n )„>q comes into play. By ([6]), 

E&* = v"P(Z' >0) 


for all n £ N. We distinguish three cases: 
CASE A. £4/1)108^(1) <0. 
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By ([TJ) and what has been pointed out in the short review of BPRE’s at the end of 
the previous section, we then infer 

l = E 4,<» = oM.^ (1) ‘ s ? 

and thus 7 < 1 , which in turn entails 

EX < E X = f < 1 

for all ?! € N as claimed. 

Case B. Eg^ (l^ogg^ (1) > 0 and ( jAsGe| >. 

Then, by (|9), 

P(z'>0) ~ c «- 3/2 (Snf E 4 (l) e ) as ?i —* 

holds true for a suitable constant c £ (0,°°) whence, by Fatou’s lemma, 


0 = limv"P(Z' >0) = liminfE^* > E (liminf X) ■ 

n—n —>00 \ n—^00 J 

Consequently, P(Surv) = 0 since inf„>o X — 1 a - s - on Sun/. 
CASEC.E4(l)log4 i (l)>0. 

Using contraposition, suppose sup n6 j - EX > 1 a ud fix an arbitrary vector a = 
(a ^)l <j<k<°° of distributions on No satisfying 

ai M) < P(x u ’ k) = x) 

for x > 1 and j,k as stated, hence 

a ii,k) > p (xU,k) = 0 ^) and £ a^ k) < P > 1 ?) 


for each ?i > 0. Possibly after enlarging the underlying probability space, we can 
then construct a BwBP (/V v , Z a v , ) v€ y coupled with and of the same kind 

as the original BwBP such that 




= ax' k] 


and 


y(j;k) ^ y{j,k) 

A a,i,v ~ A i,v 


a.s. 


for all 1 < j < k < °°, v £ V, i > 1 and x > 1. Then Z a v < Z v a.s. for all v £ V and 
since the choice of a has no affect on the cell splitting process, we have v a = V > 1 
and thus for 9 £ [0,1] 
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= E(E(Zi a |A„) 8 ) 

= £ ^ 

< E ^<E4(l)» 

1< i<k<oo v 


(13) 


where v a , Z^ XfJ’ k> and a (i have the obvious meaning. Recalling = EX^ ,k \ 
a similar calculation as in leads to 


E 1°8 8a,A q ( 1 ) < Elo gc?A 0 (!) < 0. (14) 


We are now going to specify suitable a to complete our argument. 

For M £ N let a{M) = (a44) (A4)) i<y<^ <cx , be the vector defined by 


(Xx’ k ^ (M) 



if 1 < x < M, 
if x > M , 


if Hj,k 4 1 /M, and 0Cq’^ = 1 if 11 , /. < 1/M. Then the BwBP with truncation a(M) 
satisfies the condition |AsGe]i, and we can fixMgN such that sup )ig pjE^^ M ^ n > 1, 
because ^^(m) n t as M Then, by what has already been proved under 
CASE B in combination with ( fl3| , ( [T4| and V a (M) > 1’ we infer 

oMi E *4 (1 >" £ 0 “ E *Wi)* > \ < 15 > 


which contradicts (tv). 

“(ti) => (iv)” Suppose that E£?* < 1 for all n € No and further v > 1 which, by 
< 0 >. entails lim„^ooP(Z' > 0) = 0, thus Elogg^ o (l) < 0 or Elog (1 — gAo(O)) = °°- 
We must show that E log £^(1) < 0 and info<e<i Eg^ o (l) e < v 1 . 

Assuming Elogg^l) < 0, the second condition follows from ( fTT) if ( |AsGe[ > is 
valid. Dropping the latter condition, suppose that info<e<i Eg^ o (l) e > V . Then 
we can find a M > 1 and construct a suitable “a (M) -coupling” as described above 
such that ( fl5] > holds. But since ( |AsGe| > is fulfilled for the truncated process, we arrive 
at the contradiction 

supE^* > supEtZ„* (M))! > 1 

ngN «6N 

by referring to ( fi~4| and by what has already been established for a BwBP with 
a subcritical ABPRE, i.e. Elogg^ n (l) < 0. Thus info<e<iEg^ n (l) e < v _1 holds 
even if ( |AsGe[ > fails. 

It remains to rule out that E log g' Ao (1) >0. Assuming the latter, we find with the 
help of Jensen’s inequality that 
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jnf ,l°g%Ao(l) e > n jnf0EI°g ( ^o(lj > 0 

U ^ fc/ <. 1 


or, equivalently. 


inf 

o<e<i 


E^l ) 0 



v 


(which implies info<e<i Eg^fT) 0 = 1). Use once more a suitable “a-coupling” (not 
necessarily of the previous kind) and fix a in such a way that 


1 = inf Eg' A (l) e > inf Eg' ay , (l) 0 > 
O<0<1 0 O<0<1 “ ,A) 


1 

v’ 


which implies subcriticality of the ABPRE (Z ' an ) n >0 by taking the logarithm and 
using Jensen’s inequality. As above, we thus arrive at the contradiction 

supE^* > supE.5^* n > 1 
n6N «6N 

by using the already established results for a BwBP with subcritical ABPRE. This 
completes the proof of (b). □ 


Glossary 


{Pk)k>0 

T 

T„ 

N v 

% 

%,k 

Tn 

7 —* 

n 

z v 

6¥> 

■=^n 


T* 

67 ~* 

'- y n 

xf: 


X (•>*) 

Hk 

y 

V 

(Yn)n> 0 


offspring distribution of the cell population 
cell tree in Ulam-Harris labeling 

subpopulation at time (generation) n [= {v G T : |v| = n}] 

number of daughter cells of cell v € T 

cell population size at time n [= #T„ = £| V | = „A^ V ] 

number of cells at time n containing k parasites 

the infinite vector (3T n fl, i ,...) 

the infinite vector (^, 1 ,^, 2 , ■•■) 

number of parasites in cell v 

number of parasites at time n [= £ veT Z v ] 

population of contaminated cells at time n [{v G T„ : Z v > 0}] 

number of contaminated cells at time n [= #T*] 

given that the cell v has k daughter cells v 1 vk, the j th component 

Xj J v kl of this Ng-valued random vector gives the number of offspring 

of the i th parasite in v which is shared into daughter cell v j. 

generic copy of the xj~*’ k \ i € N, v G Y with components X^' k ) 

= E xVM 

mean number of offspring per parasite [= EiZ) = Ltr>i PkYl)=\ Pj,k\ 
mean number of daughter cells per cell [= E N = Lk>iPk] 
infinite random cell line in V starting at Vo = 0 . 
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(Z')„>o associated branching process in iid random environment (ABPRE) 

A = (A„)„>o and a copy of (Z y J„> 0 . 
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